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ABSTRACT 


A  series  of  experimental  surface  pressure  tests  were  conducted  on  a 
sharp  leading  edge,  flat  plate  model.  The  tests  were  performed  in  the 
AFFOL  High  Temperature  Facility  (HTF).  Nominal  test  conditions  were 
Mach  number  10  and  a  freestream  Reynolds  number  per  foot  ranging  from 
1.65xl05  to  4.6xl05.  Angle  of  attack  ranged  from  0°  to  20®.  The  reduced 
data  are  compared  with  several  viscous  interaction  theories  and  with  test 
results  obtained  in  other  facilities. 
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SECTION  I 
INTRODUCTION 

At  hypersonic  speeds,  the  deceleration  of  the  fluid  by  the  viscous 
processes  in  a  boundary  layer  generally  produces  very  high  temperatures. 
One  result  of  these  high  temperatures  Is  an  increase  In  the  thickness  of 
the  boundary  layer  over  the  thickness  encountered  at  the  same  frees tream 
Reynolds  number  at  lower  speeds..  This  Increased  boundary  layer  thickness 
results  In  a  significant  change  In  the  effective  body  shape  which  In  turn 
generates  large  Induced  pressures.  These  pressures  are  transmitted  to 
the  Invlscld  shocklayer  thereby  altering  the  shockwave  strength  and  shapa. 
This  Interaction  between  boundary  layer  and  shockwave  Is  commonly  referred 
to  as  a  viscous  Interaction. 

Viscous  Interaction  Is  one  of  a  number  of  Important  high  speed 
phenomena  that  require  proper  evaluation  prior  to  making  any  realistic 
aerodynamic  performance  estimate  for  a  high  speed  vehicle.  It  Is  the 
purpose  of  this  report  to  experimentally  determine  the  accuracy  of 
Interaction  theories  based  on:  (1)  the  classical  hypersonic  viscous 
Interaction  parameter,  X,  originally  formulated  by  Lees  and  Probsteln 
(Reference  12),  and  (2)  the  binary  collision  scaling  parameter, 

(p«,  X*)'1^2,  due  to  the  work  of  Harney  (Reference  20). 

Data  was  obtained  by  conducting  a  series  of  wind  tunnel  tests  on  a 

sharp  leading  edge  flat  plate  model  In  the  AFFDL  High  Temperature  Facility 

(HTF).  Nominal  test  conditions  were  Mach  number  10  and  a  freestream 

5  5 

Reynolds  number  ranging  from  1.65  x  10  to  4.6  x  10  per  foot.  Angle  of 
attack  ranged  from  0°  to  20°.  The  reduced  data  are  compared  with  several 
viscous  interaction  theories  and  with  test  results  obtained  In  other 
facilities. 

The  present  study  Is  limited  to  the  classical  regime  of  viscous 
Interaction  where  continuum  theory  and  no-slip  boundary  conditions  apply. 
Neither  the  leading  edge  problem  nor  molecular  effects  are  considered. 

A  preliminary  report  on  the  test  Is  presented  in  Reference  1. 
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SECTION  II 

TEST  EQUIPMENT  DESCRIPTION 

1.  AFFDL  HIGH  TEMPERATURE  FACILITY  {HTF) 

The  HTF  Is  a  blcwdown  wind  tunnel  which  operates  at  a  nominal  Mach 
number  of  10.  Stagnaxlon  temperatures  from  2200°R  to  3500°R  are  generated 
by  an  alumina  dioxide  pehble  bed  storage  heater.  Stagnation  pressures 
may  be  varied  from  100  to  600  psla.  The  tunnel  operates  intermittently 
with  nominal  run  times  of  2  to  3  minutes*  depending  upon  the  mass  flow 
Into  the  vacuum  system.  It  has  an  open  jet  test  section  with  a  nominal 
15  Inch  diameter  test  core. 

The  basic  design  of  the  wind  tunnel  and  its  physical  characteristics 
are  presented  In  References  2  and  3.  The  nozzle  system  described  In 
these  reports  has  been  modified  and  a  new  diffuser  has  been  designed  to 
match  the  enlarged  nozzle  exit  diameter. 

The  test  was  conducted  using  a  contoured  nozzle  having  a  0. 861-Inch 
throat  diameter  and  a  24-inch  exit  diameter.  Nominal  test  section  Mach 
number  was  10.  Reservoir  pressures  were  varied  from  200  to  500  psla. 

2.  EXPERIMENTAL  FLAT  PLATE  MODEL 

A  sketch  of  the  pressure  model  showing  the  location  of  the  10  pressure 
taps  and  one  thermocouple  Is  presented  in  Figure  1. 

The  model  Is  made  of  15-4PH  stainless  steel.  The  leading  edge  was 
held  to  a  maximum  radius  of  .002  Inches.  The  top,  bottom,  and  sloping 
front  surfaces  are  1/8-inch  thick  and  are  supported  by  a  1/2-1 nch  thick 
C-frame.  Both  surfaces  are  silver  soldered  to  a  copper  water  jacket. 

The  model  was  cooled  during  the  test  to  minimize  transpiration  effects 
due  to  a  temperature  differential  between  the  model  surface  and  the 
transducer. 
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Figure  1.  Flat  Plate  Pressure  Model  Three 
Dimensional  View 
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3.  INSTRUMENTATION 

Surface  pressure?  *-^re  measured  directly  by  ten  Consolidated  Control 
Corporation  Ultrady.^t  i'r-issure  Transducers.  In  the  pressure  range  tested, 
quoted  laboratory  accmacy  Is  -  0.1%  RMS  of  the  full  range.  In  the 
present  case  this  amounts  to  t  0.0156  mm  Hg  over  the  range  0.4  to  16  mm 
Hg.  Each  transducer  was  connected  to  a  model  orifice  by  30  Inches  of 
0.140"  0.  D.  x  0.120"  I.  0.  stainless  steel  tubing  and  125  Inches  of 
0.0938  ID  copper  tubing.  Transducer  output  was  amplified  by  3  KHZ 
carrier  amplifiers  manufactured  by  Natel  Engineering  Company.  The  amp¬ 
lified  signals  were  then  digitized  by  a  Control  Data  Corporation  8032C 
analog/digital  converter  and  stored  on  magnetic  tape. 

One  Iron/constantan  thermocouple  was  used  to  measure  the  surface 
temperature.  Originally,  two  thermocouples  were  available;  however,  one 
was  damaged  after  the  model  was  assembled  and  before  It  was  tested. 
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SECTION  III 

METHOD  OF  DATA  REDUCTION 

1.  TEST  SECTION  CONDITIONS 

For  the  test  conditions  run,  it  is  permissable  to  use  the  thermally 
perfect  air  approximation  as  described  in  Reference  4.  This  method  has 
been  Incorporated  In  a  CDC  160A  computer  program.  Temperature  and 
pressure  fluctuation  corrections  to  Mach  number  have  been  estimated  from 
AEDC  wind  tunnel  correction  data  (see  References  2  and  3).  A  summary  of 
the  test  conditions  for  which  data  were  obtained  Is  presented  in  Table  I. 

2.  FLAT  PLATE  PRESSURE  MEASUREMENTS 

The  output  signals  of  the  pressure  transducer  and  thermocouple  ampli¬ 
fiers  were  digitized  and  stored  on  magnetic  tape.  A  CDC  160A  computer 
was  used  to  put  the  fundamental  data  onto  computer  cards.  The  data  were 
subsequently  used  as  input  data  to  a  CDC  6600  computer  program  that  per¬ 
formed  the  final  data  reduction  and  compared  the  experimental  results 
with  theoretical  predictions. 

3.  TEST  PROCEDURE 

After  the  wind  tunnel  reached  the  nominal  test  condition,  the  model 
was  injected  Into  the  flow  field  and  rotated  to  the  proper  angle  of 
attack.  Then  data  was  taken  at  15  sec  Intervals  until  the  tunnel  flow 
broke  down.  The  raw  data  were  examined  for  constancy  of  P  .  P..,  and 
Tw.  The  test  point  judged  to  be  most  representative  of  the  data  taken 
was  reduced  and  compared  with  the  various  viscous  interaction  theories. 

All  of  the  data  shown  In  this  report  were  taken  only  after  the  sur¬ 
face  temperature  of  the  water-cooled  model  had  equilibrated  and  the 
pressures  had  staballzed.  This  point  was  usually  reached  after  the 
model  had  been  exposed  to  the  flow  60  to  90  seconds.  Therefore,  any 
tubulation  lag  time  effects  on  the  measured  pressures  should  be  minimal. 
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However,  a  certain  amount  of  surge  In  the  wind  tunnel,  on  the  order  of 
-  U  In  PM,  did  occur.  The  effect  of  this  nonsteady  flow  Is  discussed 
In  more  detail  In  Section  V. 

The  repeatability,  consistency,  and  lack  of  scatter  In  the  data  tend 
to  substantiate  the  quoted  accuracy  of  the  pressure  transducers  and  their 
associated  amplifiers. 
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TABLE  I 


SUMMARY  OF  AFFDL  TEST  CONDITIONS 


Figure 

No, 

41nv 

(DEG) 

Po 

(PS1A) 

To 

(°R) 

Re  /FT 
x  10"® 

*** 

1-0 

3.7 

0 

197.2 

2374 

9.82 

1.68 

709 

0 

199.7 

2277 

9.90 

1.81 

669 

0 

298.7 

2422 

9.89 

2.41 

710 

0 

503.5 

2460 

9.92 

3.92 

711 

0 

504.1 

2494 

9.93 

3.81 

707 

10,11 

6.11 

496.8 

2388 

9.71 

4.34 

753 

5.11 

301.5 

2340 

9.72 

2.72 

734 

5.12 

200.9 

2334 

9.59 

1.89 

723 

12,13 

10.12 

504.4 

2331 

9.73 

4.60 

773 

10.12 

300.8 

2388 

9.73 

2.61 

793 

10.13 

201.0 

2409 

9.85 

1.65 

751 

14,15 

15.21 

297.9 

2226 

9.95 

2.78 

816 

16,17 

15.28 

502.7 

2418 

9.83 

4.15 

876*’ 

IB, 19 

11.19 

500.0 

2373 

9.88 

4.21 

770 

20,21 

12.23 

504.0 

2368 

10.04 

4.08 

828 

22,23 

13.21 

504.1 

2406 

9.81 

4.22 

797 

24,25 

14.20 

504.0 

2407 

10.00 

4.00 

870 

26,27 

17.25 

503.9 

2361 

9.86 

4.32 

852 

26,29 

19.24 

502.6 

2352 

10.01 

4.16 

851 

30,31 

20.20 

503.8 

2411 

10.02 

3.96 

862 

*  This  value  represents  the  best  estimate  for  Tw.  A  more  exact  value  Is  not 
available  because  the  thermocouple  output  saturated  the  capabilities  of  the 

CDC  8032C  analog  to  digital  converter. 
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SECTION  IV 
REVIEW  OF  THEORY 


1 .  DERIVATION  OF  DA5IC  PARAMETERS 


In  the  development  of  the  classical  weak  and  strong  Interaction 
theories i  the  need  arises  for  a  representative  value  of  py  within  the 
boundary  layer.  This  ran  be  done  once  a  suitable  reference  temperature, 
T*,  Is  established.  The  choice  of  a  suitable  reference  temperature  is 
most  Important,  and  in  the  past  many  different  methods  have  been  used  to 
determine  T*.  Chapman  and  Rubesln  (Reference  5)  proposed  using  the  wall 
temperature.  Monaghan  (Reference  6)  used  a  function  based  on  a  relation 
developed  by  Crocco  (Reference  7), 


Cheng  (Reference  8)  simplified  the  Monaghan  function  by  assuming  the 
ratio  of  the  temperature  at  the  boundary  layer  edge  to  total  temperature 
to  be  neglibible, 


T*  = 


(2) 


Eckert  (Reference  9)  chose  a  function  based  on  a  semi  empirical  derivation, 


(3) 


This  equation  reduces  exactly  to  Monaghan's  if  the  recovery  factor,  r, 
is  set  equal  to  1/1 .32. 


In  this  report  the  Eckert  reference  temperature  is  chosen  as  being 
more  meaningful  when  local  flat  plate  similarity  exists.  By  employing 
a  linear  viscosity-temperature  relation,  the  appropriate  py  ratio  is 
determined. 
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The  largest  drawback  In  applying  the  Eckert  formula  Is  the  need  to 
calculate  the  recovery  factor.  For  laminar  flow  the  recovery  factor  Is 
taken  to  be 


r*  a 


(4) 


where  the  Prandtl  number,  Pf*,  is  calculated  at  the  reference  temperature. 
Details  are  contained  In  Appendix  I  for  the  calculation  of  Pr*. 


2.  VISCOUS  INTERACTION  PRESSURE  THEORIES 


A  number  of  theories  are  available  In  the  literature  which  are 
pertinent  to  the  viscous  Interaction  problem.  The  earliest  work  Is 
contained  In  several  closely  related  studies  by  Lees  (References  10  and 
11),  Lees  and  Probsteln  (References  12  and  13),  Bertram  (Reference  14), 
and  LI  and  Nagamatsu  (Reference  15).  The  results  were  summarized  by 
Hays  and  Probsteln  (Reference  16). 

The  classical  weak  and  strong  viscous  Interaction  formulations,  based 
on  similarity  concepts  and  presented  herein,  are  from  Reference  16. 
Additional  theories,  put  forth  by  Moore  (Referen-.e  17)  and  Dewey 
(Reference  18),  based  on  the  concept  of  locally  similar  flows  are  also 
noted . 

In  addition  the  binary  collision  scaling  concept,  originally  pro¬ 
posed  by  Blrkhoff  (Reference  19)  and  applied  to  the  viscous  Interaction 
problem  by  Harney  (Reference  20)  and  by  Harney  and  Petrie  (References 
21  and  22),  Is  given.  Although  Harney  developed  the  binary  collision 
scaling  approach  to  show  that  a  meaningful  similitude  exists  in  very 
high  energy,  dissociated,  nonequlllbrlum,  frozen  air  flows,  the  theory 
Is  applicable  to  lower  energy  air  flows  and  is  presented  for  comparison. 

a.  Classical,  Second  Order,  Weak  Interaction 

t 

Tho  phenomenon  of  weak  viscous  Interaction  has  been  reported  In 
References  10  through  18.  Reference  16,  which  Is  widely  available,  gives 
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the  formulation: 

p  =  1  +  Y(^)  d°rig  J  +  Y(Y4  +  J2dorig  **  (5) 

6 

For  P  -  0.725  - 

dor18  ■  (t;)+ -145  W 

This  JQri*g  is  based  on  the  Crocco  reference  temperature.  For  an  alternate 
version  of  dQr^g  based  on  the  Eckert  reference  temperature  see  Appendix  I, 


The  (Pa/Pg)  ratio  encountered  in  Equations  7  and  8  is  to  be  determined 
from  tangent-wedge  theory  (see  Appendix  II,  Equation  II-1). 

b.  Classical,  Zeroth  Order,  Strong  Interaction  (<S^nv  =  C°) 

As  in  the  weak  interaction  case,  strong  viscous  interaction  has 
been  investigated  theoretically  in  References  10  through  18.  The  result, 
as  given  in  Reference  16,  is 

(y-  1  )Vfy(y  +  1)  G  J  .  (10) 

'  00  * 
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The  parameter  G,  which  Is  a  function  of  the  Prandtl  number,  the  specific 
heat  ratio,  and  the  temperature  ratio,  Tw/T0,  was  determined  from  the 
values  given  In  Table  9-2  of  Reference  16,  and  Is  plotted  In  Figure  2. 
These  values  are  essentially  the  same  as  those  calculated  by  LI  and 
Nagamatsu  (Reference  23).  They  are  based  on  the  assumptions  of  unit 
Prandtl  number,  y  =  1.4,  a  Falkner-Skan  pressure  gradient  function,  3, 
equal  to  0.286,  and  a  linear  viscosity  -  temperature  relation,  to  ■  1. 
Cheng  (Reference  8)  obtained  much  the  same  result  for  a  Prandtl  number 
of  3/4, 

£  *  0.332^/5’  ■  |l+2.6f^|X.  OD 

Poo 


[ + ‘-‘ft)]*  • 


Dewey  (Reference  18)  used  the  concept  of  local  similarity  to  extend 
the  range  of  viscous  Interaction  theory  to  Include  the  case  where  the 
Falkner-Skan  parameter,  e.  Is  allowed  to  vary  In  the  streamline  direction, 
and  the  Prandtl  number  is  not  necessarily  equal  to  unity.  Other  exact 
similarity  restrictions  were  relaxed  but  they  will  not  be  considered  here 
since  they  are  not  germane  tc  the  present  problem.  In  the  strong  inter¬ 
action  limit  to  Dewey's  solutions,  the  pertinent  equation  is 

£  .  3  (y  -  l)v^Y  (y  +  l)  (^I,)X  (12) 

Pco  8  ' 

where  I-|  was  obtained  from  tables  contained  In  Reference  18,  and  is  pre¬ 
sented  in  Figure  2. 

Moore  (Reference  17)  used  a  modified  local  flat  plate  similarity 
concept  to  formulate  solutions  in  both  the  weak  and  the  strong  inter¬ 
action  regions.  His  results  using  the  simplified  wall  temperature 
factor,  a,  where 


a  =  1+2.6 


(13) 


correlate  fairly  wall  with  more  exact  forms  given  in  References  16,  18, 
and  23,  and  shown  in  Figure  2.  Using  thin  relation  in  place  of  G  or 
the  strong  interaction  limit  to  the  similar  solutions  may  be 
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written  as 

£  =|  (Y-  1) 

Poo 


Vy  (y  +  i) 


1  +  2.6 


(14) 


This  corresponds  to  Dewey's  result  and  compares  favorably  with  Cheng's 
result  (Eq.  11). 


c.  Classical,  First  Order,  Strong  Interaction  (6^ny  =  0°) 

The  only  first  order  calculations  that  have  been  carried  out  are 
for  the  case  of  an  insulated  flat  plate  at  zero  degrees  angle  of  attack 
with  Pr  =■  1  and  y  “  1.4.  To  the  first  order,  the  induced  pressure  on  a 
flat  plate  is 


f  *  P(1)+  P  (o)  X  (15) 

where  and  correspond  to  the  zeroth  and  first  order  terms  in  an 
\ 

asymptotic  series  expansion  representing  the  static  pressure  distribution 
along  the  outer  edge  of  the  boundary  layer.  Calculations  for  p^  and 
p^  have  been  carried  out  by  many  authors.  The  various  results  (Refer¬ 
ence  24,  Table  I)  are  summarized  in  Table  II.  Rather  than  determine 
p^  by  solving  the  von  Karman  momentum  Integral  relation  (Reference  10), 
it  is  recognized  tnat  eventually  the  strong  interaction  region  will 
transition  to  a  weak  interaction  region  where  can  be  determined  by 
equating  the  strong  and  we«k  solutions  for  various  I  and  T./T.  ratios. 

( 1  \  WO 

For  y  =  1.4  and  Pr=  0.725,  pv  '  varies  from  about  0.90  for  ^=2  and 
Tw/Tq  =  0.2  to  about  0.98  for  J  =  0.5  and  T //TQ  =  0.8.  These  resuits 
are  comparable  to  Lees  (Reference  10)  value  of  0.92.  Rather  than  select 
any  derived  value,  is  arbitrarily  set  equal  to  1.0  since 
P/Po,  — — —  1  as  X”  — —  0.  Little  error  is  Involved  in  making  this  rather 
arbitrary  selection.  As  a  result  it  may  be  postulated  that  in  the  strong 
interaction  regime  p  /T 

£  =  1  +  3_  (y  -  1)  Vy(t  +  1)  1+2.6 

P.  16  L  Vo 


13 


AFFDL-TR-73-58 


where  the  simplified  temperature  function  Is  adopted  In  favor  of  the 
more  complex  functions  Q  and  1^. 


d.  Classical,  First  Order,  Strong  Interaction  (<5^nv>0°) 


In  the  case  of  a  sharp  flat  plate  Inclined  at  an  angle  of  attack 
with  respect  to  the  freestream,  the  strong  Interaction  result  Is 
(Reference  16) 


£ 

Poo 


P0)<M.  s1nv 

F175- 


X 


(17) 


where  p^  and  p^  are  the  same  as  for  the  zero  degree  angle  of  attack 
equations.  However,  calculations  beyond  the  zeroth  order  for  the 
functions  pv  '  and  pv  '  have  so  far  not  been  carried  out.  Therefore,  no 
comparison  between  the  first  order,  strong  Interaction  solution  and 
experimental  data  at  angles  of  attack  can  be  made. 


e.  Binary  Collision  Scaling  Interaction 

The  classical  viscous  Interaction  formulation  leads  to  the 
Probsteln-Lees  hypersonic,  viscous  Interaction  parameter,  J.  The  form¬ 
ulations  by  Blrkhoff  (Reference  19),  Harney  (Reference  20),  and  Harney 
and  Petrie  (References  21  and  22)  lead  to  a  description  of  hypersonic 
viscous  interaction  In  terms  of  a  binary  collision  scaling  parameter, 
(pjl*)’1^2.  With  this  approach  the  hypersonic,  small  disturbance 
solution  becomes  (Appendix  IT,  Equation  11-36) 


(18) 


(cp ♦5»)'  (Cp-f^y-  [I.O..H ♦r>]*  [(cp-^K  +  ^/ri 
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TABLE  II 

COMPARISON  OF  STRONG  INTERACTION  THEORY  ANALYSIS 
FOR  y  “  1.4  AND  TW/TQ  «  1 


£.  pO)  *  p(o>  y 

K  00  | 

Theory 

p<°> 

p(» 

Stewartson 

0.555 

0.489 

j  1.525 

00 

Oguchl 

0.555 

0.194  Hr.WS 

00 

Lees 

0.52 

0.92 

LI  &  Nagamatsu 

0.514 

Nagakura  &  Naruse 

0.510 

0.759 

Cheng 

Mann  &  Bradley 

0.8 

Aroesty 

-0.26  (I„3/2/M„) 

Galkin 

0.505 

-0.101  (xw3/2/mJ 

Tien 

0.528 

0.743 

Pan  &  Frobsteln 

0.651 
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Appendix  II  shows  that  the  approximate  result  (Equation  11-30)  for  the 
secondary  correlation  factor,  G^,  corresponds  very  closely  with  the  more 
complex  exact  solution,  Equation  11-28.  The  approximate  result  is 


G- 


(19) 


The  stretched  axial  coordinate  transformation,  X*,  that  produces 
local  flat  plate  similarity  In  the  presence  of  small  disturbance  pressure 
gradients  Is 

x 

X*  =  p/Idx  (20) 

o  H 

In  the  strong  Interaction  limit  X*/X  ■»  2/3,  whereas  In  the  weak  Inter¬ 
action  limit  X*  ■  X.  A  useful  linearized  approximation  for  departures 
from  the  strong  Interaction  limit  Is  (Reference  20) 


For  the  case  of  a  flat  plate,  Equation  18  reduces  to  a  quadratic 
expression  for  the  weak  and  strong  Interaction  limits,  and  to  a  cubic 
expression  for  the  linearized  approximation  for  departures  from  the 
strong  interaction  limit.  However,  as  is  pointed  out  in  Appendix  III, 
the  exact  solutions  are  not  as  satisfactory  as  the  approximate  solutions. 
Therefore,  the  approximate  solution  Is  adopted  on  the  grounds  that  it 
was  derived  from  a  consistent  set  of  assumptions. 

cp  =  [j  .015  (1  +  Y)]  1/2  G“  (PMx*)  •  1/2  (22) 


The  appropriate  X*  -  X  relation  is  applied  to  obtain  solutions  for  the 
weak,  linearized,  and  strong  interaction  limits  in  Appendix  III. 
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(a)  Weak  Interaction  Limit  Solution,  X*  ■  X 

Cp  ■  [ 1.035  (I  +  G  (/»  X)’J  (23) 

(b)  Linearized  Departure  Solution,  X*  ■  ^x[l'+^- 


Cp  -  (i5")  +  }  [l.OI5(l  +  y)]  d/ipmX)  -  (  i  )  (24) 

(c)  Strong  Interaction  Limit  Solution,  X*  *  2  X/3 

cp  •  y/\  [,  015  <l  +  r>]1  o'(^xr<  (25) 


(2)  «1nv>0« 

In  the  case  of  a  wedge  Inclined  to  the  freestream,  It  Is  not 
possible  to  discount  the  term  as  In  the  <5^nv“  0°  case  because  of  the 

-I  »2 

Importance  of  the(Cp  +  pj'q,,,  )  <S^nv  term  on  the  right  hand  side  of 
Equation  18.  Therefore,  due  to  the  complex  nature  of  this  equation,  an 
Iterative  technique  must  be  used.  This  was  accomplished  on  a  CDC  -  6600 
digital  computer. 


It  Is  apparent  from  Equations  19,  23,  24,  and  25,  that  combining  the 
wall  temperature  effect  and  binary  collision  Interaction  effect  Into  one 
unified  correlation  parameter  Is  possible.  This  parameter  Is  designated 
ff  and  Is  used  as  the  Independent  variable  In  the  plots  of  Cp  In  the 
figures. 
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SECTION  V 
DISCUSSION 

In  the  literature,  numerous  theories  for  predicting  the  effects  of 
viscous  Interaction  can  be  found.  Rather  than  attempt  to  compare  each 
theory,  many  of  which  differ  only  slightly,  this  report  limits  Itself  to 
considering  only  the  classical,  second  order,  weak  (Equation  5)  and  first 
order,  strong  (Equation  16)  Interaction  solutions  in  addition  to  the 
binary  collision  scaling  analysis.  In  both  the  classical'  and  binary 
collision  theories,  a  secondary,  temperature-dependent  correlation  para¬ 
meter  Is  encountered  In  addition  to  the  primary  correlation  parameters 
J  and  (p^X*)"^2,  respectively.  When  comparing  the  classical  analysis 
to  experiment,  the  simplified  temperature  function  1/2  a  (Equation  13)  Is 
used  In  this  report  Instead  of  the  more  complex  relations  G  and  */l  I,, 
presented  In  References  16  and  18. 

In  the  binary  collision  scaling  analysis,  a  similar  function,  desig¬ 
nated  G'  Is  encountered.  Harney  (Reference  20),  with  the  aid  of  a 
number  of  simplifying  assumptions,  was  able  to  arrive  at  a  function 
dependent  only  on  the  wall -to-total -temperature  ratio  (Equation  19).  In 
Appendix  II,  the  approximate  values  for  G'  are  <ihown  to  agree  remarkably 
well  with  the  more  exact  values  afforded  by  Equation  11-28.  Therefore, 
the  approximate  equation,  being  both  accurate  and  simple,  Is  used  In  this 
report  for  data  correlation. 

A  summary  of  the  wind  tunnel  conditions  that  were  run  in  the  AFFDL 
tests  Is  presented  In  Table  I. 

Table  III  contains  an  identification  code  to  be  used  In  relating  the 
theoretical  results  of  Section  IV  to  the  experimental  data  presented  In 
Figures  3  through  31. 

The  analysis  is  divided  Into  two  parts,  one  for  the  flat  plate  ori¬ 
ented  parallel  to  the  freestream,  and  one  for  the  flat  plate  Inclined  at 
an  angle  of  attack  to  the  freestream. 


18 


AFFDL-TR-73-58 


TABLE  III 


THEORY  IDENTIFICATION  CODE* 
FOR  USE  WITH  FIGURES  3  THROUGH  31 


Code 

Type  of  Interaction  Theory 

Equation 

Equation 

WI 

Classical,  Second  Order,  Weak 
(dorig  from  Equation  6) 

5 

5 

WIE 

Classical,  Second  Order,  Weak 
(dorig  ^rom  ERuat^on 

5 

5 

SI 

Classical,  First  Order,  Strong 

16 

- 

BCU 

Binary  Collision,  Weak 

23 

18 

BCL 

Binary  Collision,  Linearized 

24 

18 

BCS 

Binary  Collision,  Strong 

25 

18 

*The  term  "INV"  encountered  In  these  figures  refers  to  the 
inviscid,  perfect  gas  results  for  wedges  at  various  angles 
of  attack  (Reference  4). 


AFFDL-TR-7  3-58 


1.  FLAT  PLATE  PARALLEL  TO  THE  FREESTREAM  (61ny  ■  0°) 

In  Figure  3  results  of  the  AFFDL  zero  degree  angle  of  attack  tests 
are  compared  with  the  classical,  second  order,  weak  Interaction,  and 
first  order,  strong  Interaction  theories.  Kendall's  experimental  results 
(Reference  25)  are  compared  with  the  classical  theories  In  Figure  4, 
while  results  from  Bertram's  experiments  are  shown  In  Figures  5  and  6. 

In  Figures  7,  8,  and  9  the  data  of  Figures  2,4,  and  5  are  compared  with 
the  predictions  of  binary  collision  theory. 


In  the  derivation  of  the  classical  and  binary  collision  theories, 
fairly  restrictive  restraints  are  Imposed  on  the  magnitude  of  M^d.  The 
effective  angle  of  attack,  6,  Is  the  sum  of  the  geometric  angle  of 
attack,  <5.|nv,  plus  the  Induced  angle  of  attack  attributable  to  the 
boundary  layer  displacement  effect. 


6  «  6 


Inv 


.  d<$* 
+  3x 


(27) 


The  classical,  strong  pressure  Interaction  region  Is  characterized  by 
the  hypersonic,  strong  shock  limitation, 

(M^)2  »  1  (28) 


whereas  the  classical,  weak  Interaction  regime  Is  characterized  by 

(Mm  d)2  «1  (29) 

The  only  limitation  on  the  binary  collision  analysis  Is  that 


(Mwd)  ^  1 


(30) 


which  Is  less  restrictive  than  the  classical  limitations. 
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Figure  4.  Comparison  of  Classical  Interaction  Theory 
with  Experimental  Pressure  Distribution  on 
Kendall's  Flat  Plate  at  ■  0° 


Figure  5.  Comparison  of  Classical  Interaction  Theory 

with  Experimental  Pressure  Distribution  on 

Bertram's  Flat  Plate  at  *  0° 

inv 
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Figure  9.  Comparison  of  Binary  Collision  Interaction 

Theory  with  Experimental  Pressure  Coefficient 
on  Bertram's  Flat  Plate  at  b  0° 
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Similar  linearizations  In  the  binary  collision  analysis  reduce  errors 
generated  when  the  restriction  set  up  In  Equation  30  Is  violated.  The 
AFFDL  data  (Figure  3)  at  pQ  ■  500  psla  approaches  the  weak  Interaction 
limit  while  Bertram's  case  (Figure  5)  best  approximates  the  strong 
Interaction  limit.  The  rest  of  the  data  falls  between  the  two  Interaction 
regions  although  most  of  Kendall's  data  (Figure  4)  corresponds  to  the 
weak  Interaction  region.  For  different  cases  having  similar  values  of 
M*  lx  /  *  trying  degrees  of  correlation  between  data  and  theory  exist. 
In  Kendall's  case,  good  agreement  between  weak  Interaction  theory  and  the 
data  occurs  for  J  $  2  which  corresponds  to  *  0.25,2whereas 

no  similar  correlation  Is  apparent  for  Bertram  until  -  0.06, 

and  even  then  the  agreement  Is  not  as  good  as  with  Kendall.  However, 
this  Is  probably  due  to  other  reasons  which  are  discussed  further  on. 

No  particularly  good  correlation  between  AFFDL  data  and  the  weak  Inter¬ 
action  theory  exists,  although  correlation  with  the  Eckert  reference- 
temperature-dependent  theory  Is  better  than  with  the  Crocco  reference 
temperature.  In  the  region  where  weak  Interaction  effects  would  be 
expected  to  dominate  (7-  4)  the  data  Is  6X  to  7%  higher  than  the 
Eckert-based  theory,  and  7%  to  9X  higher  than  the  Crocco  based  results. 
Agreement  between  weak  Interaction  theory  and  data  is  best  when  adia¬ 
batic  wall  conditions  prevail  but  deteriorates  as  Tw/Tq-*0. 


When  the  strong  interaction  solution  (Equation  16)  Is  compared  with 
AFFDL  and  Kendall's  results  (Figures  3  and  4),  excellent  correlation 
with  experimental  results  Is  observed.  When  Bertram's  data  Is  conslderec 
the  data  trend  Is  linear  and  parallel  to  the  theory  although  a  signifi¬ 
cant  amount  below  It.  In  order  for  a  linear  pressure  correlation  to 
exist,  the  expression  for  Mw(d6*/dx)  must  be  far  more  complicated  than 
Indicated  In  Equation  32.  Since  many  approximations  and  series  trun¬ 
cations  were  made  In  the  derivation  of  Equation  32  (Reference  16),  It  Is 
quite  likely  that  In  the  Intermediate  region,  where  neither  weak  nor 
strong  Interactions  predominate,  the  expression  for  Moo(d5*/dx)  is  quite 
complex.  Bertram  had  other  "unpublished"  data  (Reference  16,  Figure 
9-4)  that  correlates  very  well  with  the  classical,  first  order,  strong 
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Interaction  theory  for  an  insulated  flat  plate.  These  data  are  repro¬ 
duced  In  Figure  6.  For  all  four  sets  of  data  It  appears  that  a  linear 
pressure  variation  does  exist  and  that  Equation  16  will  accurately  pre¬ 
dict  the  pressures  In  the  weak.  Intermediate,  and  strong  Interaction 
regions. 

In  Figure  7  the  AFFDL  data  are  compared  with  the  binary  collision 
scaling  theory.  If  the  AFFDL  data  at  pA  »  200  psla  Is  examined  In 
detail,  it  can  be  shown  that  an  115!  error  in  Mw  plus  an  additional 
3 . 3%  error  due  to  the  linear  approximation  for  Cp  occurs.  If  a  14.3% 
correction  to  the  binary  collision  theory  Is  made,  excellent  correlation 
between  data  and  strong  Interaction  theory  results.  In  contrast,  If 
similar  corrections  are  made  to  Kendall's  and  Bertram's  results,  worse 
correlation  occurs.  Similar  corrections  to  the  classical,  weak  Inter¬ 
action  theory  would  promote  better  agreement  between  data  and  theory, 
while  no  Improvement  on  classical,  strong  Interaction  theory  could 
posslbly^esult.  Although  the  ratio  between  experimental  values  for 
(mm  and  {ftp)*  determined  by  classical,  strong  Interaction 

theory  Is  relatively  constant  over  a  large  range  of  p/p^,  the  same  Is 
not  true  when  the  binary  collision  prediction  (Equation  33)  Is  ratioed 
to  experimental  values.  Clearly,  the  binary  collision  prediction  for 
(m^  for  the  AFFDL  tests  does  not  follow  the  same  trend  as  do 

either  the  classical  weak  or  strong  Interaction  predictions.  The 
deviation  from  the  trend  Indicates  that  the  binary  collision  result  for 

J  *  ^ 

3x~  9^ ves  to°  1ow  a  prediction  for  the  AFFDL  tests.  The  reason 

for  the  low  prediction  Is  due  to  the  cool  wall  temperatures  encountered 
In  the  AFFDL  tests  as  opposed  to  the  higher  temperatures  associated  with 
Kendall's  and  Bertram's  tests.  Due  to  the  Inclusion  of  a  reference 
temperature  In  conjunction  with  a  linear  velocity  profile  assumption  In 
the  determination  of  8*/6  In  Equation  IH8,  a  1%  error  In  T*/Tfi  Is 
capable  of  producing  a  6%  error  In  the  binary  collision  prediction  for 
Cp,  whereas  In  the  classical  prediction  a  15!  error  In  T*/Tg  will  cause 
less  than  a  1/25!  error  In  the  p/pw  ratio.  Reference  9  shows  that 
greater  errors  In  T*  will  occur  when  the  wall  Is  cooled  than  when  it  Is 
at  adiabatic  or  near  adiabatic  condltons.  Therefore,  the  poor  correla¬ 
tion  between  binary  collision  theory  and  AFFDL  data  is  due  to  the  strong 
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theoretical  dependence  upon  reference  temperature.  The  effect  of  a  more 
accurate  estimation  for  T*  would  show  up  only  In  Equation  11-28  and  not 
In  Equation  30,  which  Is  only  a  simplified  approximation.  Consequently, 
good  correlation  should  not  be  expected  between  current  binary  collision 
theory  and  data  obtained  at  other  than  adiabatic  wall  conditions. 

In  Figure  9  data  from  Bertram's  experiments  (Ref. 26  )  Is  compared 
with  theory.  The  data  agree  fairly  well  with  the  linearized  departure 
from  the  strong  Interaction  limit  In  the  binary  collision  scaling  theory. 
However,  the  poor  correlation  of  the  same  data  and  classical  Interaction 
theory  (Figure  5)  should  be  kept  In  mind  before  any  final  decision 
regarding  the  apparent  good  correlation  with  binary  collision  predictions 
Is  made.  If  the  source  for  this  data  (Reference  26)  Is  examined,  It  Is 
apparent  that  the  pressure  distribution  in  the  direction  of  flow  varied 
as  much  as  6-1  /25t  over  a  distance  of  3  Inches  In  the  vicinity  of  the 
nozzle  exit  plane.  It  was  In  this  region  that  the  model  test  data  were 
obtained.  In  addition,  a  pressure  difference  of  as  much  as  5%  can  be 
observed  over  a  cross  model  distance  of  only  1.2  In.  Since  Bertram's 
model  was  not  cooled  and  data  were  taken  at  less  than  adiabatic  wall 
conditions,  the  wall  temperature  was  continuously  Increasing,  so  the 
pressure  was  also  continuously  Increasing.  Therefore,  at  a  given  instant 
the  pressure  being  sensed  at  the  transducer  may  not  have  been  the  same 
as  that  being  experienced  at  the  model  surface  due  to  unfavorable 
tubulatlon  lag  time  characteristics.  For  these  reasons,  not  as  much 
creditability  should  be  conferred  on  Bertram's  data  as  on  either 
Kendall's  or  the  AFFDL  data. 

Before  going  on,  some  comment  should  be  made  with  regard  to  Kendall's 
data.  Although  the  results  were  published  In  1957,  the  pressures 
measured  were  well  within  the  Instrumentation  capabilities  at  that  time. 
Little  scatter  and  good  repeatability  may  also  be  observed  from  his  data. 
Although  Kendall's  pressure  taps  were  located  over  a  flat  plate  distance 
nearly  Identical  to  Bertram's,  no  flat  plate  side  effects  or  other 
anomalies  are  dlscernable.  The  flat  plate  used  by  Kendall  was  allowed 
to  reach  the  equilibrium  temperature  (adiabatic  wall)  before  pressure 
measurements  were  made.  Therefore,  no  pressure  lag  time  effects  would 
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be  expected.  Thus,  one  can  have  much  more  confidence  In  using  Kendall's 
data. 

2.  FLAT  PLATE  INCLINED  TO  THE  FREESTREAM  (i1nv>0°) 

For  <5.,  >0°,  weak  Interaction  effects  predominate  when  Nl  -r—  Is 
i  nv,  j  00  ox 

less  than  or  of  the  order  of  one  or  when  g^-<(5inv  with  ny 
arbitrary.  This  criteria  was  satisfied  for  all  the  AFFDL  test  cases. 

A  major  difficulty  encountered  In  evaluating  the  data  at  higher 
angles  of  attack  was  selection  of  the  most  representative  data  point. 
Although  the  HTF  Is  essentially  a  steady  state  wind  tunnel,  a  very  small 
amount  of  flow  field  fluctuation,  on  the  order  of  t  U  In  p^  occurred,, 
which  caused  the  reduced  data  to  vary  In  such  a  manner  that,  for 
51nv"  ^°0'  data  from  at  least  one  point  from  each  run  fell  on  or 
between  the  weak  and  strong  binary  collision  solutions.  It  was  this  test 
point  which  was  chosen  to  be  representative  of  all  the  data  obtained  on 
that  particular  run. 

The  measurements  of  surface  pressures  on  a  wedge  at  5,  10,  and  15 
degrees  angles  of  attack  are  shown  In  Figures  10  through  15.  It  can  be 
seen  that  the  second  order,  weak  Interaction  theory  agrees  quite  well 
with  the  experimental  surface  pressure  on  the  wedge  for  all  three  angles 
of  attack.  Agreement  between  data  and  classical  weak  Interaction  theory 
can  be  Improved  slightly  If  the  Eckert  reference  temperature  Is  used  In 
the  calculation  of  d^g.  Results  using  Eckert's  reference  temperature 
afford  approximately  a  1S6  increase  in  the  theoretical  pressure  estimate. 

At  6^nv  «  5.11°,  the  exact  solution  to  the  binary  collision  theory  suffers 
some  of  the  nonlinear  effects  mentioned  In  Appendix  III.  At  the  two 
higher  angles  (10®  and  15°),  the  weak  limit  to  the  binary  collision 
solution  provides  reasonably  good  correlation  with  data.  Moreover,  the 
classical  weak  and  binary  collision  weak  solutions  agree  remarkably 
well  with  each  other  for  the  three  angles  of  attack  shown. 


34 


AFFDL-TR-73-58 


Figure  12.  Comparison  of  Classical  Interaction  Theory  with 
Experimental  Pressure  Distribution  on  the  AFrDL 
Flat  Plate  at  *  10.12° 
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Figure  15.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  Pressure  Coefficient  on  the 
AFFDL  Flat  Plate  at  6j  =  15.21° 
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As  the  angle  of  attack  is  Increased,  the  linearized  binary  collision 
theory  increases  to  the  point  where  the  results  are  indistinguishable 
from  the  strong  Interaction  limit  values.  This  result  is  due  to  the  very 
large  value  for  p/p^  which  forces  the  second  term  on  the  right  hand  side 
of  Equation  21  toward  zero. 

The  test  run  at  pQ  ■  200  psia  and  s.|ny  *  15°  did  not  last  long 
enough  for  the  model  pressures  to  reach  equilibrium  conditions.  Conse¬ 
quently,  the  data  from  this  run  Is  not  included  In  this  report.  In 
addition  the  <5^nv  ■  16°  case  was  Inadvertently  never  run. 

Data  at  61ny  ■  15°  and  p0  *  500  psia  were  also  obtained,  but  the 
thermocouple  output  saturated  the  analog-to-digital  equipment  and  an 
accurate  wall  temperature  was  not  obtained.  However,  the  very  pronounced 
and  abrupt  departure  of  data  from  the  straight  line  prediction  at  this 
condition  Is  considered  Important  enough  to  be  Included  in  this  report 
and  is  shown  In  Figures  16  and  17. 

Originally,  it  was  planned  that  data  would  be  taken  at  pQ  ■  200,  300, 
and  500  psia,  and  <5.jny  »  5°,  10°  and  15°.  However,  the  unusual  pressure 
dip  near  the  leading  edge  at  pQ  =  500  psia  and  6^ny  a  15°  warranted  addi¬ 
tional  exploration  at  angles  of  attack  in  the  vicinity  of  15°.  These 
additional  tests  were  made  at  pQ  =  500  psia  only  due  to  limited  test 
time.  Results  for  pQ  «  500  psia  and  5^ny  =  11°,  12°,  13°,  14°,  17°,  19° 
and  20°  are  shown  In  Figures  18  through  31. 

Near  the  leading  edge  the  experimental  pressures  are  lower  than  the 
predicted  values.  The  degree  of  departure  Is  most  evident  in  the 
vicinity  of  6jny  =  15°,  when  the  freestream  unit  Reynolds  number  is 
4.15  x  105.  When  the  Reynolds  number  is  decreased,  the  effect  is  less 
pronounced.  No  departure  is  evident  at  <5jnv=  5°,  but  some  nonlinear 
effects  are  still  noticeable  at  S.  ■  20°. 
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Figure  17.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  Pressure  Coefficient  on  the 
AFFDL  Flat  Plate  at  a,..,  -  15.28° 
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Figure  18.  Comparison  of  Classical  Interaction  Theory 
with  Experimental  Pressure  Distribution  on 
the  AFFDL  Flat  Plate  at  -  11.19° 
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Figure  19.  Comparison  of  Binary  Collision  Interaction 

Theory  with  Experimental  Pressure  Coefficient 
on  the  AFFDL  Flat  Plate  at  61nv  -  11.19° 


Figure  21.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  Pressure  Coefficient  on  the 
AFFDL  Flat  Plate  at  51ny  ■  12.23° 
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Figure  25.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  Pressure  Coefficient  on  the 
AFFDL  Flat  Plate  at  «1ny  -  14.20° 
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Figure  29.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  Pressure  Coefficient  on  the 
AFFDL  Flat  Plate  at  61nv  «  19.24° 
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Figure  31.  Comparison  of  Binary  Collision  Interaction  Theory 
with  Experimental  I’ressure  Coefficient  on  the 
AFFDl  Flat  Plate  at  61ny  *  20.20° 


Other  experimenters  (Reference  27,  28,  29,  and  30)  have  attributed 
this  pressure  fall-off  near  the  leading  edge  to  merging  of  the  boundary 
layer  and  shock  wave.  In  the  present  tests,  there  Is  some  question  as 
to  whether  the  pressure  decrease  Is  due  to  merging  or  to  some  other 
unknown  phenomena,  since  the  effect  becomes  more  pronounced  with  Increas¬ 
ing  freestream  unit  Reynolds  number.  This  Is  contrary  to  the  usual  flat 
plate  explanation  based  on  the  onset  of  the  merged  layer  regime  (References 
27  and  28)  and  contrary  to  the  experimental  results  presented  In  References 
29  and  30. 

In  Reference  27  tests  on  a  sharp  flat  plate  showed  that  the  onset 
of  merging  began  at  about  ■  0.15.  However,  this  limit  Is  based  on 
the  zero  degree  angle  of  attack  situation.  Pan  and  Probsteln  (Reference 
28)  estimated  the  pressure  deviations  In  the  merged  region  with  respect 
to  a  zeroth  order,  strong  interaction  equation.  Again,  however,  the 
calculations  are  only  for  a  flat  plate  at  zero  degrees  angle  of  attack 
and  for  a  Mach  number  on  the  order  of  ten  or  greater. 

Vas,  lacavazzl,  Carlomagno  and  Bogdonoff  (Reference  29)  Investigated 
experimentally  the  effect  of  angle  of  attack  on  pressure  distributions 
on  and  near  the  merged  regime  for  sharp  leading  edge  wedges.  They  were 
not  able  to  correlate  the  merged  region  length  with  a  V  based  on  either 
freestream  or  local  Invlscld  wedge  conditions.  For  M^  ■  25,  Tw/TQ»  .15 
and  Rew  /IN  ■  9300  the  onset  of  merging  occurred  In  the  vicinity  of 

■  0.17  at  61nv  »  5°  while  at  S^nv  ■  14®  merging  began  when  .21. 

The  experiments  of  McCrosky,  Bogdonoff,  and  Genchi  (Reference  30) 
Indicate  that  the  onset  of  merging  at  angle  of  attack  is  not  a  simple 
function  of  V^,  and  that  merging  occurred  at  between  one  to  two  Inches 
for  61nv  ■  8°  and  at  about  one  Inch  for  6inv  -  15c .  This  result  is  not 
in  agreement  with  the  trend  of  the  AFFDL  results  for  PQ  *  500  psia  where 
merging  becomes  evident  at  about  one  Inch  at  v  =  10°  and  at  about 

1.2  Inch  at  6,|  »  15°.  However,  Identification  of  the  point  of  depart¬ 

ure  from  an  Interaction  type  solution  is  somewhat  arbitrary  and,  in  the 
present  case,  difficult. 
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Harney  and  Petrie  (Reference  22.  Figure  4)  encountered  merging 
effects  In  experiments  made  In  argon-air  mixtures.  Merging  effects 
became  noticeable  at  all  angles  of  attack  when  IT  reached  a  value  of 
6000  to  8000. 

Present  AFF0L  results  begin  to  deviate  from  Interaction  predictions 
for  values  of  FT  as  low  as  2000  (Figure  17)  and  7@  as  low  ?$  0.022.  No 
correlation  between  the  data  and  7  based  on  either  freestream  or  local 
Invlscld  wedge  conditions  exists.  Clearly,  the  present  results  show 
that  significant  deviations  from  the  strong  Interaction  prediction  occur 
for  values  of  FT  and  7  considerably  smaller  than  that  encountered  In 
other  similar  tests  (References  22.  29,  and  30)  for  compression  angles 
of  attack. 

At  present  no  theory  exists,  that  can  correlate  pressure  data  on  flat 
plates  at  angle  of  attack  when  merging  effects  are  present. 
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SECTION  VI 
CONCLUSIONS 

Detailed  experimental  pressures  were  generated  on  a  flat  plate  both 
at  zero  and  compression  angles  of  attack  In  a  hypersonic  flow  field. 
Pressure  data  are  compared  with  several  Interaction  theories  and  with 
data  taken  in  other  wind  tunnel  facilities. 

At  $^nv  ■  0°  the  classical  strong  Interaction  theory  agrees  with  the 
AFFDL  data  taken  under  cool  wall  conditions  much  better  than  either  the 
classical  weak  or  binary  collision  predictions.  The  poor  correlation 
between  data  and  binary  collision  theory  Is  attributable  to  the  cool 
wall  conditions  that  existed  during  the  tests. 

At  all  angles  of  attack  It  Is  shown  that  the  classical  weak  Inter¬ 
action  theory  can  be  Improved  slightly  If  Eckert's  reference  temperature 
Is  used  In  place  of  Crocco's. 

At  angle  of  attack  both  the  classical  and  binary  collision  weak 
Interaction  solutions  correlate  very  well  with  data  except  for  a  pressure 
dip  near  the  leading  edge  at  angles  of  attack  at  and  In  the  vicinity  of 
15°.  It  Is  thought  that  this  dip  Is  due  to  the  merging  of  the  boundary 
layer,  although  values  for  the  rarefaction  parameter  do  not  satisfy  the 
usual  conditions  necessary  for  the  onset  of  merging. 
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APPENDIX  I 

DETAILS  OF  CALCULATIONS  FOR  THE  PRANDTL  NUMBER,  REFERENCE 
TEMPERATURE,  AND  CHAPMAN-RUBESIN  VISCOSITY-TEMPERATURE  RELATION 
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1.  PRANDTL  NUMBER,  Pr 

A  least  square  curve  fit  of  data  presented  In  Reference  32  1 
to  determine  the  Prandtl  number, 

Pr  -  A  +  U  T  +  C  T2  , 

where 

200aR  i  T  S  1000°R 
A  -  >8238997 
B  --.2975599  X  10-3 
C,  •  .1522722  X  10-6 
and 

1000°R<T  i  1800oR 
A  ■  ,6772683 
B  --.145707  X  10“4 
C  -  .1579909  X  10’7 

2.  REFERENCE  TEMPERATURE,  T* 

Following  Eckert  (Reference  9)  the  relation  for  T*  Is 


T*  ■  %  ♦  -  T.)  *  .22  (Tm  -  Te) 

To  find  the  adiabatic  wall  temperature 

T  ■  T  +  r*  (T  -  T  ) 
aw  e  '  o  e' 

It  Is  necessary  to  evaluate  the  recovery  factor,  r*,  which  for  a 
laminar  boundary  layer  Is 
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Finally,  an  expression  for  T*  Is  obtained 

T*  -  .5(Te  +  Tw)  +  .22  (T0  -  Te)  (1-5) 

Since  the  ?r  -  T  approximation  Is  limited  to  a  second  order  ciJrve 
fit.  It  Is  possible  to  solve  for  T*  explicitly. 


For  1000°R  <  T*  i  1800°R 

A  -  1  -  7.6467498  (TQ  -  Tfi)2  X  10“8 
B  ■  7.0522188  (T0  -  Te)Z  X  10"7  -  (Tfi  +  T^) 

C  -  .25  (Tfi  +  Tw)2  -  3.2779785  (TQ  -  Te)2  X  10  ‘2 
For  200°R  <  T*  *  10008R 

A  -  1  -  7.3699744  (Tq  -  Te)2  X  10’9 

B  ■  1.4401899  {Tq  -  Tfi)2  X  10'5  -  (Te  +  Tj 

C  -  .25  (T#  +  Tw)  -  3.9876745  (TQ  -  Tg)2  X  10  'Z 

3.  CHAPMAN  -  RUBESIN  RATIO,  C* 

When  T*  is  known,  It  Is  a  simple  matter  to  calculate  the  reference 
viscosity,  y*  by  using  either  the  Hlrschfelder  (Reference  33  )  or 
Sutherland  Equations. 

For  T*  <200°R 

u  *  -  8.051  T*  X  10"10  ( -L-B-  \ 


(1-7) 
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For  T*  >  200°R 

w*  "  2,27  Tr'^~T?gTF  X  10"8  (  (I-8) 


By  Incorporating  the  usual  approximation  that  the  pressure  Is  constant 
across  the  boundary  layer,  the  Chapman-Rubesin  ratio, 


c*  .  fiinl 

pewe  * 

(1-9) 

becomes 

v* 

c*  “  * 

(1-10) 

4.  ECKERT  REFERENCE  TEMPERATURE  RELATION  FOR  dQr1g 

The  Eckert  reference  temperature  can  be  used  in  conjunction  with 
the  zero  pressure  gradient  result  for  boundary  layer  displacement  thick¬ 
ness  growth, 


to  obtain 


S*  *  1  .7208x/  ^yRe* 


6* 


p  U 
ie 


1.7208 


(1-11) 


(1-12) 


After  differentiating  this  expression  to  find  d<5*/  dx,  Equation  1-2  Is 
used  to  get 


where 


^orig 


.4302 


d6*  .  J  . 

3T  “  jr  dor1g 

v 


(1-13) 


1  + 


r 


—  +  (  Y-  1)  (.44  yPr*  + 


.  (1-14) 
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APPENDIX  II 

A  SUMMARY  OF  THE  HYPERSONIC  VISCOUS  INTERACTION  PHENOMENA  IN 
TERMS  OF  HARNEY'S  BINARY  COLLISION  SCALING  PARAMETER,  (p„x*)“ 
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The  analysis  contained  In  this  appendix  Is  entirely  attributable  to 
Dr.  D.  J.  Harney  of  the  AFFDL  and  Is  presented  here  for  the  convenience 
of  the  reader.  Because  It  circumvents  the  need  to  have  an  accurate  free- 
stream  Mach  number,  as  Is  required  In  classical  theory.  It  Is  most  valu¬ 
able  In  connection  with  experiments  In  very  high  energy  wind  tunnels 
where  non-equilibrium  effects  become  significant  and  Mach  number  deter¬ 
mination  Is  uncertain.  It  is  included  in  this  report  because  there  Is 
no  good  reason  why  the  theory  should  not  be  applicable  to  tests  conducted 
In  equilibrium  flow  fields.  The  method  can  be  simply  stated  as  fo1 1  '"'s : 
The  pressure  coefficient,  Cp,  can  be  represented  as  a  function  of  the 
gas  properties,  the  frsestream  flow  conditions,  and  the  effective  angle 
of  attack  of  the  body.  The  effective  angle  of  attack  Is  the  sum  of  the 
geometric  angle  of  attack  plus  the  added  body  thickness  caused  by  the 
displacement  effect  of  the  viscous  boundary  layer  as  It  exists  in  the 
vUcous  Interaction  region.  The  hypersonic,  small  disturbance,  tangent- 
wedge  theory  due  to  Llnnell  (Reference  34 )  Is  used  to  account  for  the 
Invlscld  processes,  while  the  compressible,  laminar  boundary  layer  effects 
are  accounted  for  by  a  reference  property  transformation  of  Blaslus* 
Incompressible  flat  plate  boundary  layer  equations.  The  pressure  grad¬ 
ient  resulting  from  the  Interaction  effect  is  accounted  for  by  a  suitable 
transformation  of  the  streamwise  coordinate.  This  analysis,  attributable 
to  Harney  (Reference  20)  Is  described  here. 

1.  INVISCID  EFFECTS 

For  planar  flow  the  small  disturbance,  tangent-wedge  theory  of 
Llnnell  (Reference  35)  Is 


where  In  our  case  Sis  the  sum  of  the  geometric  angle  of  attack  plus  the 
viscous  Induced  angle  of  attack  attributable  to  boundary  layer  displace¬ 
ment  effects.  ns* 

5  ■  «1nv  *  M'  ■  <n'2! 
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Within  the  realm  of  small  disturbance  theury  it  is  also  assumed  that 
M  Sind  »  1 ,  so  that 

OQ 


K  ( l±L\\  (  2  * 

.\8 .  r+|  +  i  /  2  \t/(Tr+ 1\ 

y  v  a  ) + 

1/  2  ^  l  M.Sin  8  //  V  2  / 

(n-3) 


By  substituting  the  first  two  terms  of  the  binomial  expansion  given  above 
into  the  small  disturbance  equation,  the  following  result  (Reference  31, 
Equation  5)  is  arrived  at: 

Sin*8  <y  +  U  +  ty  +  l)  M„*  S»n*8  '  (II-4) 

Using  the  usual  perfect  gas  relationship  for  M^, 

*2  f*m  “•  "  i  •  (II-5) 


the  pressure  coefficient  expression  reduces  to 


CP 

Sin*  8 


ly  ♦  •> 


q.  Sin1 8 


(H-6) 


Because  the  second  term  of  this  expression  is  much  less  than  the  first, 
the  effect  of  y  is  contained  primarily  in  the  first  term. 

ine  coefficient 


y+l 


(11-7) 


may  then  be  considered  a  constant,  b. 

Cp  bp_  v 

~TZ  ■  (y-H)  + - ~ 

5m  8  q  S I  n* £■ 


or,  rewriting 


p  -(l  +  bl1  p„ 

— 


(y  +  »)  Sin* 8 


(H-8) 


(11-9) 


67 


AFFDL-TR-73-58 


Harney  (Reference  31)  defined  the  left  hand  side  of  this  equation  as  C' 


Cp  *  (yfll  Sin1 8 


(11-10) 


or,  In  terms  of  C  , 

P 


CD  •  c0  + 


•  cB  + 


(11-11) 


The  best  justification  for  this  expression  is  that  it  approximates  the 
exact  and  Linnell's  results  for  MMS1n6  *2  quite  well  with  b  ■  1.  If  b 
Is  allowed  to  vary  so  that  a  better  fit  to  exact  results  is  obtained, 
then  b  ■  3/4  (Reference  3l).  and 


(y-H)  Sin* 8 


■  I.OI5 


(11-12) 


These  "constants"  are  entirely  empirical. 

For  a  wedge,  the  inviscid  pressure  coefficient  is  given  as 


Cp  •  1.015  (y  +  l)  Sin*  8  J-  ^ 


(H-13) 


2.  COMPRESSIBILITY  EFFECTS 


Using  a  reference  temperature  concept,  the  Blasius  result  for  an 
Incompressible  boundary  layer  may  be  applied  to  the  compressible  bound¬ 
ary  layer  problem,  thereby  making  it  possible  to  calculate  the  local 
slope  of  the  boundary  layer.  If  and  only  if  similarity  exists  can  It 
be  shown  that 

(„_i4) 

dx  e  At  e  2 


where  8  Is  the  momentum  thickness  and  Cf  is  the  local  coefficient 
of  friction  e 


C*.  "  n 


(11-15) 
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Employing  the  classical  Blaslus  formulation  for  Cf  and  incorporating 

e 

the  reference  property  concept,  Harney  (Reference  20)  gets,  for  the  zero 
pressure  gradient  case, 

•  0.332  ^  (R*tx)  ‘  (IM6) 

where 

R**x  *  /»*«**//£  (II_17) 


as 


The  displacement  thickness,  momentum  thickness  ratio  can  be  expressed 


(11-18) 


or  by  using  a  reference  property  concept  In  conjunction  with  a  linear 
velocity  profile 


where 


(II-19) 


(11-20) 


This  may  be  compared  with  the  Blasius  result  for  (6*/e)^nc.  which  is 
2.591.  Harney  (Reference  20,  Equation  A-8)  chose  to  use  a  value  of 
2.61  rather  than  either  the  linear  or  the  Blasius  values.  It  is  shown 
in  part  4  that  by  using  the  2.61  factor  in  the  calculation  of  the 
approximate  temperature  function,  GC  a  good  correlation  results  with  the 
exact  linear  velocity  profile  analysis  which  uses  the  3.00  value. 
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Combining  Equations  11-14,  11-16,  11-19,  and  11-20,  It  can  be  shown 

that 

4  ’  1992 ['-(&)] ("*••«)'*  (II-2,) 

3.  PRESSURE  INTERACTION  EFFECT 

Harney  (Reference  20)  accounts  for  the  nonzero  pressure  gradient 
effect,  arising  from  the  Interaction  of  the  viscous  and  Inviscld  flow 
fields,  by  a  stretching  transformation  of  the  streamwlse  coordinate. 


o 


This  transformation  produces  local,  flat  plate  similarity  In  the  presence 
of  small  disturbance  pressure  gradients.  It  Is  this  stretched  coordinate 
which  Is  used  In  place  of  the  geometric  coordinate  in  the  reference 
Reynolds  number 

#  *  P*  ut  ** 

R#«x  *  R,««*  “  "T#—  (n-23> 

By  using  the  zeroth  order  weak  and  strong  Interaction  solutions  as  limits 
on  the  pressure  variations,  Harney  found  that  X*  -  2X/3  In  the  strong 
limit  and  X*  ■  X  in  the  weak  limit.  For  departures  from  the  strong 
Interaction  limit,  he  used  a  linearized  departure  approximation 
(Reference  20,  Equation  37) 


70 


AFFDL-TR-73-58 


£ 

i 

I 


4.  INVISCID,  COMPRESSIBLE  AND  INFRACTION  EFFECTS 

Combining  the  results  from  parts  1,  2,  and  3,  It  can  be  shown  that 


(11-25) 


and,  after  some  manipulation,  that 


^  {ir +  ( ■ ^ )  <°“ '*>  W  (  Jt)]} 

(11-26) 

Here  the  Eckert  reference  temperature  (Reference  9)  and  Sutherland 
viscosity-temperature  relations  have  been  used.  If  this  expression  is 
combined  with  Equations  11-21  and  11-25, then 


#  •  ‘UP '*«*’•* 


(11-27) 


where 


3.046X10- r.  /T.\j  |('+Tp  ^/y-l\r.l8,j„/M| 

A  t  198.6  t  V2T»/J|  Mf«  V  '  2  71  VT0/J 


( i 1-28) 


Here,  the  usual  recovery  factor  approximation 

r*  ■  V? 

has  been  made. 


(11-29) 
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Harney,  by  making  a  number  of  simplifying  assumptions,  arrived  at  the 
expression* 


8 

G  -  3.84  X  I0'4  [  1  +  2.6  (  I*- )]*  (11-30) 

The  more  exact  result  for  G^Equatlon  11-28)  Indicates  that  It  Is  depend¬ 
ent  upon  far  more  than  just  the  T./T  ratio.  To  Investigate  this  point 
further,  a  plot  of  the  exact  versus  the  approximate  G'  factor  Is  shown 
In  Figure  32.  Considering  Kendall's  test  (Reference  25),  the  approximate 
solution  differs  from  the  exact  solution  by  only  -0.94%,  while  the  maxi¬ 
mum  G'devlatlon  for  Bertram's  test  (Reference  25)  Is  3.25%,  Tho  approxi¬ 
mate  value  for  G'based  on  the  AFFDL  data  differs  by  a  maximum  of  4.14% 
from  the  exact  solution,  although  for  the  data  shown  In  Figure  7  the 
deviation  Is  less  than  1%.  .A  typical  case  with  TQ  «  2000°F  and  Mm  ■  10 
Is  also  shown  In  Figure  32.  The  maximum  deviation  from  the  exact  solu¬ 
tion  Is  7%  at  TW/TQ  -  0.8. 

The  approximate  G"  relation  (Equation  11-30) agrees  remarkably  well 
with  the  more  exact  solution  (Equation  11-28)  when  the  Mach  number  Is 
low  and  the  TW/TQ  ratio  approaches  the  adiabatic  wall  limit,  and  when 
the  Mach  number  Is  high  and  the  17T  ratio  Is  small.  When  the  Mach 

W  0 

number  Is  high  and  the  T./T.  ratio  approaches  the  adiabatic  wall  limit 
the  correlation  between  the  approximate  and  the  exact  values  deteriorates 
somewhat  but  In  general  remains  quite  accurate.  Regardless,  the  maximum 
deviation  of  the  approximate  Cp  from  the  exact  Cp  Is  less  than  0.7%  for 
all  the  AFFDL  data.  Consequently,  in  this  report  the  approximate 
theory,  being  much  simpler,  is  used  for  data  correlation. 


♦Harney's  analysis  also  included  the  effects  of  a  nonequilibrium  free- 
stream.  Since  these  effects  are  not  encountered  in  the  experiments 
reported  herein,  they  arc  not  included  in  this  report. 
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Equation  11-27  may  now  be  combined  with  Equation  11-13  to  get 


Cp  ■  ••OI5(l  +  y)(8|f,v  +  ^),+  |.5= 

rk 


or 


E 


r  5  Pj* 
CP  T  q» 


1.019(1  +  f) 


Sinv  +  e' (^) 


Rearranging, 


(H-31) 


(11-32) 


CP“ 


T 


f  •  [l.OIBll+y|]^  [( 


Cp+ 


(11-33) 


For  small  angles  of  attack,  the  small  disturbance  assumption  can  be 
Invoked 

£  “  I  (H-34) 


while  for  moderate  angles  of  attack  the  right  hand  side  of  Equation  II- 
33  Is  dominated  by  the 


(11-35) 


term,  and  the  effect  of  the  small  error  Introduced  by  assuming 
u  /u^  ■  1  is  significant.  With  this  approximation,  Equation  11-33 
can  be  written  as 


(cp+5;)l(cp-|-^),-[l.0,5(,+  y)]i[(cp+5j)i8inv  +  G(/»|>*4‘r*]  . 
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5.  FLAT  PLATE  ORIENTED  PARALLEL  TO  THE  FREESTREAM  («1f)V  «  O'3) 

For  the  zero  angle  of  attack  case,  Equation  11-36  simplifies  even 
further 

(CB+Ss)i  (cp-l-lf  )*  (11-37) 


This  equation  reduces  to  a  quadratic  expression  for  the  weak  and  strong 
Interaction  limits,  and  to  a  cubic  expression  for  the  linearized  approxi¬ 
mation  for  departures  from  the  strong  Interaction  limit.  A  more  detailed 
analysis  of  the  validity  of  the  exact  solution  In  comparison  with  an 
approximate  solution  Is  given  In  Appendix  III,  where  It  Is  shown  that  the 
simplified,  approximate  equation 

t  , 

Cp  «  [l.OI5<l  +  y)],cfl/>(|Jt*f*  (11-38) 


gives  a  more  reasonable  result  than  the  exact  solution. 

6.  FLAT  PLATE  INCLINED  TO  THE  FREESTREAM  (51ny>  0°) 

For  the  nonzero  angle  of  attack  case,  Equation  11-36  has  to  be 
solved  by  Iteration. 


75 


APPENDIX  III 


A  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTIONS  FOR  THE 
BINARY  SCALING  PRESSURE  COEFFICIENT  ON  SHARP  FLAT  PLATES 
AT  ZERO  DEGREES  ANGLE  OF  ATTACK 
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APPENDIX  III 


A  COMPARISON  OF  EXACT  AND  APPROXIMATE  SOLUTIONS  Fn»  thf 
HNARV  SCALING  PRESSURE  COEFFICIENT  ON  ffi°Fuf°PufIs 
AT  ZERO  DEGREES  ANGLE  OF  ATTACK 
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1 .  EXACT  SOLUTIONS 

The  most  general  form  of  the  expression  for  the  pressure  coefficient 
on  sharp,  flat  plates  at  zero  degrees  angle  of  attack  Is  (Appendix  II, 
Equation  11-36) 


which  may  be  squared  to  give 

+  ~ r>] ■  o  (iii-z) 

This  equation  can  be  solved  for  Cp  by  employing  the  quadratic  expression 


iVfm)  •<* 


1 1  -  3 ) 


For  the  weak  interaction  limit  X*  -  X,  and  the  strong  interaction 
limit,  X*  a  2X/3,  the  solutions  are  straightforward;  however,  for  the 
linearized  approximation  about  the  strong  interaction  limit 

»*•  4 »[•+*•(*?£)] . 
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A  cubic  solution  Is  necessary.  By  substituting  the  linearized  expression 
for  X*  Into  Equation  1 1 1-3  and  rearranging,  we  arrive  at  the  cubic 
equation 


cp+ 2(iJf)  c!~ 


(II1-6) 


It  can  be  shown  that  this  equation  has  three  real  and  unequal  roots. 


The  correct  root  Is 


where 


cp  •  2  v/~^'Co‘  H)  ‘“i’c 


(in-6) 


( II I  - 7 ) 


b  «  < 


[l.OISd  *fr  y)J  G* 


ipmt) 


-(&)(&)' 


{ill) 


( n  i-8) 


(IU-9) 


and 


Co*  4> 


-y/  ytifj 


(III-10) 
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2.  APPROXIMATE  SOLUTIONS 

Employing  the  basic  hypersonic  assumption,  Mm»1,  the  left  hand  side 
of  Equation  III-l  becomes 

(Cp  +  lf)^  (cp~  Otf  “  cp  (in-n) 

The  validity  of  this  approximation  can  be  ascertained  from  Figure  33. 
With  this  approximation  Equation  III-l  simplifies  considerably 

Cp  »  [l.OI»(l  +  y)]*c'(/aj«*ri  (III-12) 

For  the  weak  and  strong  Interaction  limits  the  solutions  again  are 
straightforward: 

a.  Weak  Interaction  limit,  X*  -  X 

cp  •  (i n-13) 

b.  Strong  interaction  limit,  X*  -  2X/3 

Cp  '^|)[».OI5(l  +  yl]^G(^xri  (111-14) 

For  the  linearized  departure  from  the  strong  interaction  limit,  a 
simple  quadratic  equation  results: 


79 


AFFDL-TR-73-58 


3.  EVALUATION  OF  THE  EXACT  SOLUTION 

A  comparison  of  the  exact  and  approximate  solutions  are  shown  in 
Figure  34  for  three  nearly  Identical  wind  tunnel  conditions.  On  this 
figure  is  superimposed  a  line  representing  data  sets  from  three  separate 
wind  tunnel  tests.  The  data  and  the  approximate  solutions  tend  to 
Cp  =  0  as  the  scaling  factor,  (pjO-^2,  becomes  very  small.  The 
approximate  solution  provides  a  better  correlation  with  the  data  and 
with  what  one  would  expect  a  priori  simply  because  the  assumptions 
Involved  In  the  derivation  of  the  approximate  equation  are  consistent 
with  the  assumptions  Involved  In  obtaining  the  solution  to  this  equation. 
As  the  Mach  number  becomes  larger  and  larger  the  exact  solution  will 
approach  the  basic  hypersonic,  straight  line  result.  (See  also  Figure 
33).  Only  the  approximate  straight  line  solution  will  be  considered  In 
this  report. 
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